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We use the mean-field approximation to simplify the master equation 
for sympathetic cooling of Bosons. For the mean single-particle occupation 
numbers, this approach yields the same equations as the factorization assump- 
tion introduced in an erlier paper. The stationary or equilibrium solution of 
the resulting master equation for the one-body density matrix shows that 
the mean-field approximation breaks down whenever the fraction of conden- 
sate Bosons exceeds ten percent or so of the total. Using group-theoretical 
methods, we also solve the time-dependent master equation for the one-body 
density matrix. Given the time dependence of the mean single-particle occu- 
pation numbers, this solution is obtained by quadratures. It tends asymptot- 
ically towards the equilibrium solution. 



I. INTRODUCTION 

Sympathetic cooling of Bosons or Fermions is an important tool in Atom Optics. It is 
used in cases where the interaction between atoms is too weak for evaporative cooling to 
work. This fact calls for the development of a theory describing the process. A master 
equation for sympathetic cooling was formulated in Ref. This equation turned out to be 
too unwieldy for practical purposes. A significant simplification was introduced in Refs. . 
Two independent approximation schemes - microcanonical averaging, and a factorization 
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assumption - both led to equations that could be solved numerically. The results were in 
good agreement with each other [Q. 

In this paper, we introduce and analyze another approximation to the master equation 
of Ref . . This is the mean-field approximation often used in many-body theory. We show 
that this approach leads to the same rate equation for the mean occupation numbers as 
the factorization assumption introduced in Refs. [@]. In this sense, the mean-field approach 
contains and generalizes the factorization assumption. It allows us to study not only the time 
evolution of the mean occupation numbers but that of the full one-body density matrix. 
We do so by using a group-theoretical technique in von-Neumann space. We are able 
to construct analytically both the equilibrium and the non-equilibrium solutions of the 
master equation for the one-body density matrix. This leads to a complete theoretical 
understanding of equilibrium and non-equilibrium properties of the solutions of the master 
equation and, at the same time, yields insight into the limitations of the mean-field approach 
itself. We focus attention on the case of Bosons which is theoretically more complex and, 
therefore, more interesting. 

In Section |T| we formulate the master equation and use decoherence to reduce it to di- 
agonal form. The mean-field approximation is introduced in Section |TI|. Our algebraic 



technique is introduced in Section We then investigate the equilibrium solution (Sec- 
tion 0). While this investigation could have been done without using the algebra developed 
in Section our method shows its full power when we diagonalize the rate operator in 



Section VI and when we construct the time-dependent solutions of the master equation in 



Section VII. Section VIII contains a brief summary. 



II. THE MASTER EQUATION FOR THE DENSITY MATRIX AND 

DECOHERENCE 



We follow the notation of Refs. |1],|2|. System A is subject to sympathetic cooling and 
consists of A^^ Bosons. The master equation for the time-dependence of the reduced density 
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matrix pA{t) has the form 



dpAjt) 
dt 



-I Ha + H[ 



(1) 



Here, Ha is the sum of the single-particle Hamiltonians for the atoms in system A containing 
the kinetic energy and the harmonic-oscillator potential of the trap, while H'j^^j^ contains 
the interaction between the atoms in system A. The action of the Liouvillean C on the 
reduced density matrix pA(t) is given by 



it to say that they account for the interaction between particles in system A and those in 
the cooling system B. The latter is assumed to be at thermal equilibrium at all times. 
The labels m and n refer to the three quantum numbers {mx,my,mz) defining the single- 
particle eigenstates of an isotropic harmonic oscillator in three dimensions representing the 
trap potential, with associated creation and annihilation operators at and a^, respectively. 
Eqs. (0) and (0) are derived from the von Neumann equation for the density matrix of the 
combined system A plus B under the following assumptions: The process is Markovian, 
the correlation time for the interaction between systems A and B is much shorter than the 
cooling time, and a rotating-wave approximation applies. Equations of the form of Eqs. (|1],0) 
are often used to discuss Boson condensation phenomena 

It was shown in Ref. [Q] that if the number Nb of particles in the bath is very large, 
decoherence acts very quickly and reduces the density matrix to diagonal form. Then, Eq. 
takes the form 




(2) 



The rate coefficients F^^T are given in Ref. |1| ; we do not repeat the definition here. Suffice 



CpA = E r^;S(2ajn 



afiPA(t)aia,n - alarnal^afiPAit) 




(3) 
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III. MEAN-FIELD APPROXIMATION 



Eq. d^) serves as the starting point for our mean-field approximation. We apply this 
approximation in standard fashion by replacing on the right-hand side of Eq. @ one pair 
of creation and annihilation operators referring to the same single-particle state m or n by 
its expectation value. We use the notation 

(ata^) = tr(ata^pA) = A^m • (4) 

The quantities A'^ are the mean occupation numbers of the orbital m. For the term CpA, 
this procedure yields 



+r^"(^n + 1) (2arf,pA{t)al - alarnpA{t) - pA{t)arnal^ ^ ^^)) • 
Eqs.([l|,|) and (|) imply 

dt ^ ' 

In what follows, we explore the algebraic structure of Eq. We show that this equation 
is linear in the generators of the group SU(1,1). In the algebraic theory of non-autonomous 
quantum systems it is shown that any such linear dynamical system is integrable and, 
thus, can be solved analytically. 

IV. ALGEBRAIC STRUCTURE OF THE MEAN FIELD EQUATION 

In the first two Subsections, we introduce the mathematical tools which are used in 
Subsection [IV C| to explore the structure of Eq. (|^). 



A. Right and Left Algebras 

Let \N)ffi stand for the A^-fold occupied single-particle state m. In a Hilbert-space basis 
defined in terms of such states, the creation and annihilation operators at and a^fi appearing 



in Eq. (|^) may operate from the left on the ket state \N)rn or from the right on the bra state 
rfi{N\. In the latter case, we attach an upper index r to the operators and write them to 
the left of the density matrix pA- We do so with the understanding that for a product of m 
operators Oi x O2 x . . . x Om, we have 

0[ X X . . . X O'^pa = pAOm X . . . X O2 X Oi . (7) 

This convention (the reversal of order as we bring the factors from one side of pA to the 
other) is the cause for the differences in signs between the first and the second set of Eqs. (|TU|) 
below. As a result, the pair of operators at and generates both the left and and the 
right representations of the group hw(4). We define accordingly 

if the operators act from the left on the ket state |A^),fi; and 

if the operators act from the right on the bra state rn{N\. Together with the unit operator, 
the three operators appearing on the left-hand sides of Eqs. (||) (of Eqs. @) constitute the 
left representation hw(4)/ (the right representation hw(4)j. of hw(4), respectively). These 
operators obey the following commutation relations. 



hw{4)i : 










- -at 






= +1 . 


hw{4:)r : 




= -(^rn 










«^"^] 


= -1 . 



It is evident that hw(4)i (hw(4)r) is isomorphic (antiisomorphic, respectively) to hw(4). 
Since hw(4); and hw(4)j, act in different spaces, the operators belonging to these two different 
groups commute with each other. We write this somewhat symbolically as 
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[hw{A)i, hw{A\] = 



(11) 



B. Composite Algebra 

In order to explore the algebraic structure of Eq. (^), it is not sufficient to list the basic 
algebras hw(4)i and hw(4)r. It is neccessary, in addition, to display the structure of the 
composite algebra C obtained by joining these two algebras. The composite algebra C 
consists of the elements 

C : {nj^,n^, aJi^a^, a^a^} . (12) 
The elements of C obey the commutation relations 

["'ml ^m^m\ "'rn^'m ' 

r I't r ''til r I 

[^:n,^!n]=0. (13) 

These relations follow from Eqs. ([10|) and (|ll|). The algebra C is the direct sum 
U(1)©SU(1,1) of two commuting parts, a radical U(l) and a simple Lie algebra SU(1,1). 
Here 

Uil):{n^ = n'^-n:^} (14) 

and 

SU{1, 1) : = lin^^ + nl), = a%a^^, Kz^ = a^-^a'^} . (15) 
The elements of SU(1,1) obey the relations 
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[K^,Kt] = 2Kl. (16) 
We list the action of the four elements of C on the projectors 11^ = |A^)^^(A^|. 



KzYil = iVn^-i . (17) 

Since p{t) is diagonal in energy representation, we are concerned exclusively with the pro- 
jectors n^, i.e., with the eigenvalue (—1) of the operator ra^. Our algebraic method is not 
confined to this case, however. This is shown in the Appendix. 

C. The Structure of Eq. (||) 

We define 

J- m J- n,m^^ri i 

r"^ = E ^%{Nn + 1) (18) 

and use Eqs. ( [T^JT3| ) and ([TTp. The master equation takes the form of a rate equation, 

= Y.{2T^K^ + 2T^KZ - 2(r^ + T^)Kl - (F^ - T^)]p^{t) . (19) 

Eq. (|l^) shows that aside from the mean values A^ contained in the rate coefficients F"^ 
and r^, the time-dependence of pA{t) is given by a sum of terms each referring to one of 
the orbital subspaces only. Moreover, each such term is linear in the generators of SU(1,1) 
and is, thus, integrable. This makes it possible to solve Eq. (|19|) algebraically. 

We reduce the rate equation Eq. (|I^) by taking a partial trace. We use a Hilbert-space 
representation in terms of products of states jA'rra)- We take the trace of pA(t) over all states 



n with n ^ m and denote the result by Prfiit). From Eq. ([T9|) , we find 



= [2T^K^ + 2T^KZ - 2(r^ + T^)K'. - {T^ - T^)]prn{t) . (20) 

The omission of the interaction term HIa-a from Eq. (^) was justified in Refs. [^. We 
recall that pAif) is diagonal in energy representation. It follows that can be written as 
a sum of the projectors 11^, 

Na 

pm{t) = E ^^n^ • (21) 

Af=0 

We use this form in Eq. ( 120|) and find for < Na 

r\P^(f) 

= 2Tr,,NP^-' + 2r™(iv + - 2(r^ + r"^)ivp^ - 2r^p^ . (22) 

Eq. (p2D is similar to the time-dependent Hartree-Fock equation used in Nuclear Physics [0 . 
From Eqs. (^,^) it follows that the mean particle number A^^(t) in orbit m obeys the 



equation 

dAr^(t) 
dt 



2V^{N^ + 1) - 2r'^iV^ . (23) 



Using the definitions for the F's in Eqs. (18) we see that Eq. (^) coincides with Eq. (39) 
of the second of Refs. [0. In particular, we have dY.m^m/(^'t = so that particle number 
A'^ = ^rh is conserved. This shows that the factorization assumption used in Refs. 
is equivalent to the mean-field approximation. 

V. THERMODYNAMIC EQUILIBRIUM: STATIONARY SOLUTION 

In thermodynamic equilibrium, the reduced density matrix becomes independent of time, 
and the stationary solution of the rate equation Eq. ( ^2]) obeys dpA{t)/dt = or 

Y.{2T^K± + 2T^KZ - 2{Trn + T'^)Kl - (F^ - F™}pA = . (24) 

fh 

Inserting the expansion Eq. (^), we find the following set of equations. 

2F^iVP^-i + 2F™(iV + l)P^+i - 2(F^ + F'^)iVP^ - 2F^P^ = if iV < iV^ ; 
2F^iV^P^^-i - 2(F^ + F")iV^P^^ - 2F^P^^ = ; 

2F^(iV^ + 1)P^^ = . (25) 
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It is easily seen that these equations only have the trivial solution = for all < A^^. 
This fact shows that a consistent non-trivial solution of the mean-field equations does not 
exist if we impose on pA the constraint expressed by Eq. (^). Therefore, we modify this 
equation by writing 

oo 

Prnit) = E ^In^ • (26) 

7V=0 

All we now require is that the sum ZI^^tv^ negligibly small. Then the first of Eqs. (^5]) 

applies for all N. The solution is 

= where x = ^ • (27) 

From Eq. (^) we see that in the stationary case we always have x < 1- The normalization 
condition yields = (1 — x). The mean value A''^ is given by A''^ = — x). Conversely, 
we may replace x everywhere by Nrh/{Nffi + 1). To discuss the validity of the mean-field 
solution, we impose the constraint that P^'^/P^ < exp(— a). For Na ^ 1, this yields 
N^/Na < l/o. This shows that for the condensate (the only case where we expect A'^ 
to take sizeable values), the mean-field approximation (as defined in the framework of this 
paper) begins to fail whenever the ratio N^a/Na grows beyond ten percent or so. 

The P^ ~ x^ form a sequence which decreases monotonically with increasing A^. This 
behavior is in marked contrast to that of the equilibrium solution found first by Scully 
and reproduced, in the present context, in Ref. 0|. Scully's solution is not based upon the 
mean-field approximation. We believe that it describes correctly the equilibrium behavior 
of the fully developed condensate. 

VI. DIAGONALIZATION OF THE RATE OPERATOR 

We turn to the time dependence of the solutions of the rate equation. As a preparatory 
step, we determine in this Section the eigenvalues and eigenvectors of the rate operator F, 
defined as the operator which appears on the right-hand side of Eq. (|2^), 
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r = [2T,nKt + 2T^Kz - 2{T^ + T^)K^^ - (F^ - T^)] . (28) 
For simplicity, we omit reference to the single-particle state m, define 

a = 2TrA, b = 2r™ (29) 

and have 

T = aKt + bKZ-{a + b)K'^-^ia-b) . (30) 
We note that F is not self-adjoint, and that 

Ft = aKz + bK+ - (a + b)K'^ - ^(a - 6) • (31) 



The mean-field equation ( pOD takes the form 



Tprn . (32) 



dt 

The equilibrium solution determined in section obviously corresponds to the eigenvalue 
zero of the rate operator. The eigenfunction is given by Eqs. (|2^j27| ). We now construct the 
complete set of eigenvalues and eigenvectors of the rate operator. A similar procedure has 
also been used in Quantum Optics 0. 

To this end, it is useful to first look at a simpler problem. Given the three components 
Ji, J2, J3 of spin which obey the cyclic commutation relations [Ji, J2] = iJs, we ask for the 
eigenvalues and eigenvectors of the operator 

= aJ2 + bJ3 . (33) 

The constants a, b are real. Writing a = ^JaF^^^ sin0 and b = ^/ + b"^ cos 9, we have 



O = Va?+¥ (sin ^ J2 + cos J3) , (34) 

and we recognize immediately that O can be diagonalized by a rotation in the (2, 3) plane 
by the angle —9, i.e., by the unitary transformation U = exp{—i9Ji), 
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O = UOU^ = Va^ + Js ■ (35) 

This result can be verified algebraically. The eigenvalues of O are proportional to those of 
the operator J3, the eigenvectors are obtained by applying the unitary transformation f/^ to 
the eigenvectors of J3. The diagonalizing matrix U is not unique: Replacing 6 hj 6 ±tt also 
yields a diagonal form, with J3 replaced by (— J3). 

The form of the rate operator in Eq. (|30|) is somewhat more complicated than that of the 
operator O. Moreover, the underlying group is not SU(2) but SU(1,1). Therefore, it takes 
a similarity rather than a unitary transformation to diagonalize F. Geometrically speaking, 
the vector F is not confined to the (2, 3) plane. This makes it necessary to define a pair of 
similarity transformations, 

T_ = exp(— a.i^'") , T+ = exp(— a+i^'^) . (36) 

Eqs. ([T6| ) imply the identities 

exp{-a_K-)K+ exp(a_K") = - 2a. - a^_K~ ; 

exp(-a_ir-)ir°exp(a„is:-) = - a^K' ; 

exp{-a+K+)K- exp(a+ir+) = + 2a+K^ + alK+ ; 

exp(-a+K+) exp(a+K+) = + a+K+ . (37) 

It is easy to see that the product 

T = (38) 

diagonalizes F (in the sense that TTT'^ is diagonal) provided the coefficients a_ and a+ 
obey the equations 

a + 6 a+ - (a + 6) a+ = ; 6 (1 - 2a+a_) + (a + 6) a_ = . (39) 
The transformed rate operator reads 

TTT-^ = [2ba+ - (a + b)] K^-^{a-b) . (40) 

11 



The eigenvectors of TrT~^ belonging to the eigenvalue —1 of the operator n^n, are those 
of and have the form 11", with n = 0,1, 2,.... The right-hand eigenvectors of the rate 
operator F are correspondingly given by c„T~^ 11" with c„ a normalization constant. 

Eqs. (^) have two solutions, yielding two sets of eigenvalues 7^ for the rate operator. 
This fact is due to the same type of ambiguity as found for the operator O in the last but 
one paragraph. With n = 0, 1, 2, . . ., the solutions are 

(a) a+ = x , a- = , In = n{a - h) ; 

(6) a+ = 1 , a- = , 7n = (n + 1)(6 - a) . (41) 

For the right-hand eigenvectors of F, this yields 

Pn = CnT-' n" = Cn exp{a+K+) exp{a.K_) H" , (42) 

where c„ is a normalization constant. Eq. (|4^ ) shows that the solutions (b) of Eq. (|4l|) are 
not admissible. Indeed, we have exp{a^K^) 11° = Y.n=o a+n^, and a corresponding result 
when the operator exp(a;+-ft'"'') acts on 11^ with > 1. For solution (a) of Eq. (^Tj), where 
a+ = X < 1) the resulting series behaves as required while for solution (6) where = 1, 
the coefficients multiplying FI^ do not decrease with decreasing and, thus, violate the 
condition formulated in Section |V| for the viability of the mean-field solution. Solution (a) 
does contain the zero eigenvalue discussed in Section |V|. Moreover, all non-zero eigenvalues 
7„ of solution (a) are negative, as is to be expected on physical grounds. 

Eq. (PH) shows that the operator TFT~^ is self-adjoint and, thus, equal to [T^^yV'^T^ . 
This implies that F^ has the same eigenvalues as F and is diagonalized by the similarity 
transformation (T^^Y . The right-hand eigenvectors of F"!" are given by c„rt n". They 
coincide with the left-hand eigenvectors p„ of F. It is easy to show that the right- and 
left-hand eigenvectors of F form a biorthogonal set , 

T^iPnPm) = 5nm ■ (43) 
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VII. TIME DEPENDENT SOLUTIONS OF THE MEAN FIELD EQUATIONS 



In the fully time-dependent mean-field equations Eq. (p2D, the rate operator T depends 
upon time. This is because F is defined in terms of the coefficients a and b which in turn 
depend on time via the mean occupation numbers A^^, see Eqs.(p9D and ([T8|). We now 
show that given the solutions N^{t) of the rate equations Eq. and, thus, the time 
dependence of the coefficients a{t) and 6(t), the solutions to the full equations are 
obtained by quadratures. We also show that these solutions tend asymptotically towards 
the equilibrium solution found in Section |V|. As in Section |VJ, we suppress the label for 
the single-particle state m. We retain the symbol x introduced in Section to denote the 
asymptotic or equilibrium value of the ratio a/h. 

We assume that at time t = 0, we are given the initial density matrix 



oo 



p(0) = ^ pV (44) 

Af=0 



with the proviso that for > A^^, the coefficients are negligible. To solve Eq. (|3^), we 
define the time-dependent similarity transformation 



p[t) = T{t) p{t) = exp[-a_(t) K-] exp[-a+(t) K+] p{t) . (45) 
The time-dependent coefficients a±(t) are solutions of the differential equations 

p ' = a + - (a + 6) «+ ; ^ ' = 6 (1 - 2a+a_) + (a + 6) «_ . (46) 

Lit/ Lit' 



Comparison with Eqs. ( pStpOD shows that T(t) provides the time-dependent generalization 



of the similarity transformation T used in Section N\ As initial conditions for the differential 



equations (|4q) , we choose a±(0) = or p(0) = p(0). Using the identities Eq. ([37|), we see 



that the transformed density matrix pit) obeys 



^ = [7(t)i^°-^(a-fe)]p(t), (47) 



where 
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7(t) = 2ba+ - (a + 6) . (48) 

Integrating this equation and using the relation inverse to Eq. (^), we obtain 

p{t) = exp[a+{t)K^] exp[a4t)K-] exp(|J dr[7(r);^'° - ^(a(r) - 6(r))]) p(0) . (49) 

We have accomphshed our first goal and shown that given the time dependence of the 
coefficients a and b, the time dependence of the density matrix is obtained by quadratures. 

We now show that p(t) as given by Eq. ( ^9|) tends asymptotically (t oo) to the 
equilibrium solution established in Section 0. To this end, we first investigate the asymptotic 
behavior of a;±(t). The first of Eqs. (^61) reads da+{t)/dt = b{a^ — a/b){a+ ~ !)• We recall 
that a > 0,6 > and distinguish two cases, (a) We have a/b < 1 for all times t. Then, 
da+(t)/dt > 0(< 0) if < a+(t) < a/b (if a/b < a+(t) < 1, respectively). With the initial 
condition a+ = 0, we see that a+ approaches the value a{(X))/b{oo) = % < 1 asymptotically 
from below. This implies 7(t) ^ (a — 6) < as t ^ cxd. (b) We have a/b > 1 for some time 
interval, although a/6 ^ x < 1 must hold asymptotically. As long as a/b > 1, the initial 
condition a_|_ = implies a+(t) < 1. As soon as a/b intersects the value unity from above, 
Q:+{t) has a negative (positive) slope if it is larger (smaller) than a/b. In either case, a+ again 
tends toward x < 1 asymptotically, and 7 attains the asymptotic value a(oo) — 6(00) < as 
in case (a). To study the asymptotic behavior of a_, we define f{t) = a_ exp[/Q dr7(r)]. 
The time derivative of f{t) is given by feexp[/Q dr7(r)]. Since b is bounded and 7 negative 
for large t, this derivative tends to zero and, hence, / towards a constant /q. This implies 
that a_ diverges asymptotically. 



We turn to the asymptotic behavior of p{t). We use Eqs. (^) and (|T7D to write 



p{t) = exp[a+{t)K+] exp[(l/2) /*dr(7(r) - (a(r) - 6(r))] 

Jo 

00 ^ /N\ 

xE/W'^^'^WE («-)-^n^- (50) 

N=0 p=0 \P / 

With the help of the asymptotic behavior of f(t),a+, and a_ established in the previous 
paragraph, this becomes asymptotically 
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Pit) explxK-^] exp[(l/2) / dr(7(r) - (a(r) - 6(r))](E /(O)^P^(O)) • (51) 

As 7 tends asymptotically towards (a(oo) — 6(00)), the integral in the exponent exists. 
Actually we ought to show that the difference 7 — (a — 6) vanishes asymptotically at least 
as strongly as t~^. This we have not done. We know, however, from Eq. ( PBD that both 
tr p and the mean value iV^ of the number operator remain finite for t ^ 00. Therefore, 
the normalization constant of our solution cannot diverge asymptotically. Together with 
{En=o foP^iO)), the exponential can be lumped into a new constant C. This yields 

p{t)^Cexp[xK+]U' . (52) 

This shows that asymptotically, p becomes proportional to the equilibrium solution found 
in Section |V|. 

VIII. SUMMARY AND CONCLUSIONS 

We have used the mean-field approximation to simplify the master equation for sym- 
pathetic cooling of Bosons. We have shown that the factorization assumption of Refs. 
is equivalent to this approximation. Studying the stationary or equilibrium solution of the 
resulting master equation, we have shown that the mean-field approximation to sympathetic 
cooling is expected to break down whenever the fraction of condensate Bosons exceeds ten 
percent or so of the total Boson number in the cooled gas. This conclusion is supported by the 
observation that the equilibrium solution for the one-body density matrix differs markedly 
from the form derived for the same quantity for the case of a fully developed condensate. 
Using group-theoretical methods, we have solved the time-dependent master equation for 
the one-body density matrix. Given the time-dependence of the mean single-particle occu- 
pation numbers, the solution is obtained by quadratures. It tends asymptotically [t 00) 
towards the equilibrium solution. For the time-dependence of the fully developed conden- 
sate, quantum fluctuations are important. These can probably be calculated along the lines 
described in Section V of the second of Refs. 0. 
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We expect that the mean-field approximation will play an important role in studies of 
sympathetic cooling. In this paper, we believe to have given a comprehensive theoretical 
treatment of this approximation. 
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APPENDIX A: GENERAL APPROACH TO EQUILIBRIUM 

We remove the restriction to the eigenvalues —1 of the operator and define 

= \nUrn{k\ . (Al) 

Eqs. ([T7| ) are replaced by 

rirnlQ^ = {n-k- l)H;i^ , 



K-^T^t = nr^''^-^ . (A2) 

The eigenvalues jn,k of the transformed operator in Eq. ( ^OD are given by 

-fn,k = ^{a - b){n + k + 1) - ^{a - b) , (A3) 

where n and k run over non-negative integers and the eigenvectors are c^^fcH"''^. The right- 
hand eigenvectors of F are 

Pn,k = c„,feT-^H"''= . (A4) 
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The approach to equihbrium can also be generahzed by considering 

oo oo 

pW = EE^"''(o)n^''W (A5) 

n=Ofc=0 

where P"'*^ is defined by the initial condition. In analogy with Eq. (0), Eq. ( |A5| ) can be 
rewriten as 

pit) = exp[a+i^+] exp[(i) rfr(7(r) - (a(r) - 6(r))] exp[f dtW)] 

We have assumed here k < n for simplicity. Given the asymptotic limit of the time- 
dependent function in the equation above one can immediately conclude that as t oo, 
only the term corresponding to n = k = will survive, yielding the equilibrium state. 
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